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ABSTRACT 


Shipboard measurements of temperature and velocity fluctuations were 
performed to determine optical propagation properties of the marine 
boundary layer. Empirical expressions describing the temperature struc- 
ture parameter, Be in terms of the Richardson number, Ri, overland were 
used to evaluate data obtained for open ocean conditions. Profiles of 
the dissipation of turbulent kinetic energy, ¢«, with respect to height, 
Z, and with respect to a stability parameter, Z/L, were examined for 
Open ocean conditions. 

In general, there was little correlation between the measured c, 
and the stability parameter. However, the distribution of e with height 


for both the stable and unstable cases showed little deviation from that 


expected for near neutral conditions. 
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I. INTRODUCTION 


The advent of more complex and complete numerical models for analyz- 
ing and predicting large scale atmospheric motion has also led to better 
Specifications of the boundary layer. The boundary layer is that region 
in the first kilometer over the sea which is defined principally by 
turbulent transfer of momentum and heat. 

Optical propagation through the atmosphere is affected by the re- 
fractive nature of the medium. In addition to the regular variation of 
atmospheric refractive index with height, there exist small inhomogenei- 
ties in the refractive index associated with fluctuations in the tempera- 
ture and velocity of the air. These fluctuations cause random phase 
and amplitude distortions in propagating wave fronts and thus degrade 
Spatial and temporal coherence in the transmission. The magnitude of 
the effects places limitations on optical system performance and must be 
included in design and execution considerations. 

Descriptions of the small scale fluctuations which affect optical 
propagation have not been as complete nor in the quantity for the over- 
water regime as for the overland regime. 

Initial experimental efforts to verify turbulence theory predictions 
were conducted over land. Measurements of wind speed, temperature, and 
humidity in those investigations were made from a stable platform with 
auxiliary instrumentation well protected from external influences. 

Under these conditions relatively accurate measurements of wind speed, 
temperature and humidity were performed under a variety of weather condi- 


tions. 


Ll 





The marine environment on the other hand presents a number of prob- 
lems, especially in accessibility, platform stability, and sensor pro- 
tection. Overwater descriptions are necessary, even though considerable 
progress has been made in overland investigations. The necessity exists 
because of the increasing evidence of the influence on atmospheric 
motions by features unique to the oceanic environment. 

It is the purpose of this paper to compare overland and overwater 
results, also, to evaluate the assumption that near neutral conditions 
apply for both stable and unstable conditions when evaluating the varia- 


tion of turbulent kinetic energy dissipation with height. 
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TI. THEORETICAL BACKGROUND 


A. GENERAL 

Turbulence properties of interest in optical propagation are those 
which could be used to describe the intensity of fluctuations in the 
refractive index and that which could be used to describe the size or 
Scale of the refractive index inhomogeneities. The former has been re- 
lated empirically to scintillation, beam wander and beam spread, and 
the latter along with the former to the image resolution which is im- 
portant with respect to seeing conditions. Similarity expressions re- 
lating mean conditions to properties of the turbulent regimes and hence 
the temperature structure function and the dissipation rate of turbulent 
kinetic energy have evolved through a series of hypotheses and asso- 
Ciated measurements. 

An important consideration within turbulence theory is whether the 
turbulence is anisotropic or isotropic. In general, anisotropy exists 
within large scale eddies and isotropy exists within small scale eddies. 
The universal formulae considered will be valid only for isotropic 
turbulence. 

Anisotropic turbulence is very empirical. Every time the boundary 
changes, the turbulence changes. However at high wave numbers, which 
corresponds to small scale eddies, the turbulence should be independent 
of the boundaries. Kolmogorov (1941) postulated that at high wave 
numbers turbulence should be isotropic even though it is generated 
from and embedded in anisotropic eddies of low wave numbers. If 


Reynolds numbers are high enough, turbulence will adjust through inertial 
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transfer and viscous dissipation until a statistical equilibrium, inde- 
pendent of initial conditions, is attained. This equilibrium consists 
of energy transfer into the high wave number range that is equal to 
energy lost through viscous dissipation. 

On the basis of the isotropic nature of small scale fluctuations, 
Only one parameter is necessary to describe the intensity of the refrac- 
tive index fluctuation over many scales. It is the refractive index 
Structure function parameter, a defined as 


.* = [<n(x) - n(xtr)>J2/r°/ > fal) 


where n(x) and n(xtr) are refractive indices at two points ona 
line oriented normal to the mean wind direction and separated by the 
distance r which is less than the outer scale, Lo? the lower end 
of the inertial subrange, and greater than the inner scale, ae the 
Smallest scale of naturally occurring turbulence. 

A parallel expression, which defines the temperature structure 
function parameter is 
Ceo = ET (x) = Ti (xtr)>]yr°/> 
where T'(x) and T'(xtr) are temperature fluctuations at two 
points separated by the distance r. The refractive index is deter- 
minted primarily by density fluctuations and can, therefore, be 
related to temperature fluctuations, neglecting himidity. 


-* is related to Cs by the following 


= = [79.x107© p/t?7? ox (a 
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Fortunately, both ae and ae are readily measurable by optical 


and meteorological means, respectively.:n alternate relationship for 


2 
- 


lent kinetic energy, €, and temperature variance, x, is 


c » Which involves measurement of the rates of dissipation of turbu- 


oan = 
C. = Bye 


vs (4) 
where 8 iS an empirical constant with a value of 3.2. This last form 
enables indirect estimates of oe to be made from mean conditions 
Since e€ and yx are easily related to boundary layer fluxes and, 
eventually, profiles if steady horizontally homogeneous conditions exist. 

Expressions which relate ay to mean properties of the boundary 
layer are desirable because the small scale measurements are impractical 
to obtain in most operational or tactical regimes and situations. 

Since turbulence is nearly synonymous with temperature fluctuations, 
it is ultimately desirable to describe mean thermal stratification in 
terms of atmospheric bulk stability parameters such as the Monin-Obukhov 
length (L). In this regard, measurements of both atmospheric mean 
profiles and fluxes are required for a complete determination of atmos- 
pheric transmission behavior. Profile (36/92 and 93U/az) and bound- 


ary flux (U, and T,) parameters appear in the following expressions for 


Ri and L 
996/92 
Ri = ————= (5a) 
6(9U/3z) 
_ a, 2 
L = 6U, / KT 4, (5 /o)) 





The following similarity predictions for the dependence of ec and 
x on momentum and heat fluxes and height were considered by Wyngaard, 
et al. (1971) in deriving an empirical expression for estimating ae 


from mean stability parameters 


e2/U,” = f,(Z/L) (6a) 
pa hn UR F,(Z/L) (6b) 


The form of the empirical expression for i 1s obtained by direct 
substitution of equations (6a and b) into equation (4) yielding 


Oe Ty 8, (Z/L) (7) 


ay 
Furthermore since (Z/L) and Ri can be functionally related, 
Businger, et al. (1971), a parallel dependence on Ri can be obtained, 
Walz. 


2 _ 74/3 


C (38/32) “$4 (Ri) (8) 


+ 

The functions f : f, ‘ f. , and 3 in equations (6), (7), and 
(8) are empirical and formulated on the basis of observations of both 
temperature fluctuations, momentum and heat fluxes, and mean gradients 
of wind speed, temperature, and humidity. 

The relation expressed by equation (8) provides a desired depen- 
dence of ie on more readily measured mean stability (Z, 36/3z, and 
Ri). The forms of fr. and b3 and the data obtained from the exten- 


Sive AFCRL study of turbulence structure over a flat, unobstructed 
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Kansas plain are presented in Wyngaard. As will be shown later, avail- 
able marine data does not appear to agree for fr, as well as expected 
with the overland predictions, i.e. Hughes (1976). The overland pre- 


dictions for f, and f, in equation (6), were as follows 


4.9(1-7(Z/L))7°/3, Z/L < 0 


f,(Z/L) 


F,(Z/L) Poe OC Zl) ) 5 Z/L > 0 (9) 


Smal] scale velocity fluctuation properties are of interest in 
optical propagation because image resolution has been empirically re- 


lated to the innerscale, Ros which is defined 


= (y3ye)'/" (10) 
where y iS the kinematic molecular viscosity and e is the dissipa- 
tion rate of turbulent kinetic energy similar to c,° . € can be ob- 
tained from either one-dimensional velocity variance spectral estimates 
in the inertial subrange or from velocity structure function estimates. 
€ and he can be functionally related to mean profile and flux 
estimates (U,, T, and Ri or Z/L) on the basis of the empirical expres- 
Sion presented previously (equation 6a). For the purpose of examining 
overwater ce results from different stability conditions described by 


Richardson numbers (Ri), equation (6a) can be expressed as follows 
So ; 
e = U,~/Kz o, (Ri) Gap} 


Since Z/L and Ri are functionally related. 
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Overwater e values can also be evaluated on the basis of U, 
values estimated from them. Such estimates can be evaluated by compar- 
ing them with U, estimates computed from mean wind profile measure- 
ments. Both U, estimates utilize the following expression for the 


mean wind gradient 
aU/3z = U,/kzZ $, (Ri) iz) 


where o, (Ri) is defined in equation (11). o, (Ri) is equal to 1.0 
under neutral conditions. U, is assumed to be a constant with height 
in this expression so it is restricted to the constant surface layer. 


Integration of equation (12) for near neutral conditions yields 
U= U,/k an Z/Z, (13) 


where UB 1S assumed equal to zero, Zo is the roughness parameter. 
The logarithmic profile suggested by equation (13) is altered by 
Stability, shown in Figure 1. Zo can be eliminated from equation (13) 


by selecting mean winds at two appropriate levels so that 
U, = k(U, - U,)/2n(Z,/Z,) (14) 


It is important to note that this expression relating U, to mean wind 
Values (U, and U,) at two levels (Z, and Z,) is applicable only for 
near neutral conditions (Ri= 0). 

In near neutral conditions, turbulent kinetic energy production is 


assumed to be equal to the rate of molecular dissipation of turbulent 
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kinetic energy and the following relation is 


valid 
a 2 
e= WU, (e0/az) (15) 


combining equations (12) and (15), assuming neutral conditions >, (Ri) 


= 1.0, and solving for U, , yields 
u, = (ekz)!/3 (16) 


Thus, under neutral conditions, the friction velocity (U,) can be esti- 
mated from either mean wind profiles using equation (14) or from 
fluctuation data (involving turbulent energy dissipation) using equa- 
tion (16). 

From determined values of U, , it is also possible to calculate a 
momentum drag coefficient, Ci» corresponding to a given height in the 


Surface layer 
c= U,°/0 Cr?) 


Numerous studies have been conducted to determine a representative 
value of C, for 10 meters, Cardone (1969). Check calculations of 
Cio employing U, determined from both mean wind profiles and dissi- 


pation rates are performed in this study. 


B. STABILITY CONSIDERATIONS 
Observational experiments by Businger et al. (1971) yielded a 


definite relationship between the Richardson number, Ri, equation (5a), 
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and the Monin-Obukhov length, L, equation (5b). Figure 2 from Businger, 
et al. (1971) lillustrates this relationship quite well. The following 
expressions are approximations of the relationship between Z/L and Ri, 
proposed by Dyer and Hicks (1970) and Webb (1970) for unstable and 


Stable conditions, respectively, 


Z/L = Ri (18) 
Z/L = Ri/1-akRi (19) 


where a iS an empirically derived constant equal to 5.0. 

It is important to note that Businger, et al. (1971) observed the 
Richardson number to approach a critical value of 0.21 as Z/L 
approached 2 Suggesting that as stability increases, the flow 
becomes essentially non-turbulent. This implies that for the surface 
layer the effect of mechanical turbulence becomes negligible. 

The effect of hydrostatic stability on surface layer mechanical 
turbulence is further illustrated in Figure 3 from Businger, et al. 
(1971). Recalling equation (12) in light of Figure 3, it becomes 
apparent that the function o,(Z/L) vaires relatively little with 
respect to Z/L for unstable conditions approaching a value of 0.5. 
However, the rapid increase of o(Z/L) as the atmosphere becomes 
hydrostatically stable indicates that the near neutral assumptions 


are invoked under only slightly stable conditions. 
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A. PLATFORM AND LOCATION 

All observations were made aboard the R/V Acania while anchored in 
about 30 fathoms of water off Pt. Pinos in Monterey Bay. Open ocean 
differs from land in the effects of wave action on turbulence, in the 
nature of aerosols and fog, and humidity fluctuations. These conditions 
could best be obtained far at sea; however the cost of such activities 
made it desirable to work near land and the optical propagation experi- 
ments could be coincident. Pt. Pinos and other locations along the 
Monterey Bay shoreline provided very nearly the ideal situations for 


both the optical and meteorological experiments, Figure (4). 


B. INSTRUMENTATION AND PROCESSING EQUIPMENT 

Measurements were made at four levels. Figure 5 shows positions 
of the towers which support the instrumentation, length of cable runs, 
and position of processing and recording equipment. 

The mean wind measurements were obtained with Thornwaite Associates 
cup anemometer wind profile register systems, model number 104. Three 
cups, plastic cones reinforced with aluminum frames, are mounted on 
Stainless steel tubes attached to the main shaft at 120° intervals. When 
the cups rotate due to the wind a slotted shaft in the anemometer serves 
as a shutter between a light source and a photoelectric cell. The three 
cup anemometers have the characteristics of low starting speeds with a 
small amount of internal friction which aids in checking inertial over- 


shoot. 
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Optical ranges and R/V Acania position (normal) in Monterey 
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Velocity fluctuation measurements were obtained using TSI model 1210 
probes with sliding support shields and tungsten wires. The shields per- 
mitted isolation of the sensing area for the determination of the undis- 
turbed velocity, Vo: voltage reading before and after each experiment. 
The platinum coated tungstenwire was small enough to resolve the viscous 
dissipation scale without making corrections for wire length. The wire 
had plated ends for isolating the sensing area and thereby minimizing 
flow disturbance. Figure (6) 

Electronics associated with each probe were a TSI model 10548 
linearized anemometer and a TSI model 1056 variable decade module. The 
anemometer had a linear frequency response from (dead calm) O KHz to 
10 KHz and the variable decade module operated with a 0-60 ohm range. 

Temperature fluctuations were measured using similar size sensors. 

A platinum wire was used in place of a tungsten wire. The temperature 
fluctuation system was designed for a resolution of 0.001°C at frequen- 
cies up to 1 KHz. 

The temperature fluctuations were measured using a bridge developed 
by GTE Sylvania, the GTE Sylvania model 140. The system was slightly 
modified for use aboard the R/V Acania. 

The baseband portion of this system is basically a balanced wheat- 
Stone bridge excited by a 3 KHz signal with a synchronous detector on 
the output. Segments of a very small diameter platinum wire serve as 
temperature sensors in opposite arms of the bridge. The resistance 
temperature coefficients result in an output from the bridge which is 
proportional to the temperature difference between the two probes. The 


+ 


sensor wire is 0.5 centimeters long and 2.5x10 ° centimeters in diameter. 


This extremely small mass allows a response to temperature variations 
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Figure 6. 
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of up to 1 KHz, while electronic amplification allows temperature differ- 
ences as small as 0.004°C to be observed. 

Both wind and temperature fluctuation data were recorded on a San- 
born model 3950 fourteen channel tape recorder. Real time readout on 
an eight channel chart recorder, brush model 240, was used to check the 
quality of the signals coming from the sensors. The chart data was also 
utilized in temperature variance analyses.. 

Hewlett Packard model HP-23850 temperature sensitive quartz crystal 
probes were used to obtain mean temperature at the sea surface and four 
tower levels. Figure (7) The RF signal from the probes and a reference 
oscillator were mixed in a Hewlett Packard HP-2801A readout unit to pro- 
duce a beat frequency whose signature can be analyzed to within 0.001 
°C/Hz. Each sensor received pre-experiment calibration against a 
platinum resistance wire thermometer in a temperature controlled circulat- 
ing water bath. It was checked over the range of expected temperatures. 
The accuracy in achieving a 0.005°C correction factor was a constant 
for each probe. The tower mounted sensors were housed in aspirated 
shelters to eliminate radiation effects. Figure (8) 

Mean relative humidity information was obtained using Dunmore type 
lithium chloride sensors. Figure (9) This sensor was also placed in 
the aspirated shelter. The basic principle of operation of this sensor 
is resistance change in an electrolytic solution generating a reference 
voltage variation proportional to relative humidity change. Automatic 
temperature compensation in the instruments meet the following specifi- 
cations for relative humidity, + 3% humidity below 90% and + 4% humidity 
above 90%. Sensor calibration was accomplished by a comparative method 


uSing a Saturated saline solution in a closed container. 
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Figure 9. Aspirated shelter. Figure 10. Dunmore type 
lithium chloride humidity 


sensor. 
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Figure 8. Hewlett Packard model HP-2850 temperature sensitive 
quartz crystal probes. 
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Sensor placement required exceptionally long cable runs. Here de- 
crease in system frequency response due to their length had little 
effect in the frequency band of interest. Adjustments were made in 
the bridges for resistance and capacitance of the wire length. 

Data logging during the experiments was accomplished using the NPS 
developed MIDAS (Microprogrammable Integrated Data Acquisition System). 
This micro-processor based data acquisition system utilizes an Intel 
8008 central processor to control the sampling, averaging, and record- 
ing of mean meteorological data. All software is written in PL/M to 
facilitate the writing of a self-documenting program. 

The operator is interfaced with the system via teletype for full 
duplex input/output communications and program control over the sample 
Start time and the number of samples to be averaged before outputting. 
The operator may also alter the present preset sample list by adding 
or deleting various sensors as they come on line or become inoperable 
or inconsistent. Once initiated, the system is fully automated to 
sample the tailored list of sensors every 30 seconds and periodically 
print output values averaged over the selected interval of from one 
minute to one hour. In this study a 10 minute averaging interval was 
selected. 

Output values were printed out on the teletype in a columnized for- 
mat with the time of print as a leader. The teletype has a paper punch 
incorporated which may be activated by the operator to produce a data 
copy concurrent with the printout. This paper punch was also utilized 
to produce data cards for use with the IBM 360 digital computer. A 


magnetic cassette tape recorder has been integrated into the system 
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as a third data output device. This cassette can be interfaced to an 
HP 9830 portable computer so that profile and gradient flux estimates 
can be performed on board automatically using the basic programming 


capability of the HP 9830. 
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IV. ANALYSES PROCEDURES 


A. ANALYSES OF MEAN DATA 

Mean wind data for ten minute intervals were available from the 
Naval Postgraduate School developed MIDAS (Microprogrammable Integrated 
Data Acquisition System) output. MIDAS read values at ten second in- 
tervals and averaged them for a ten minute period. Analyses were per- 
formed on data recorded on teletype printout and also on the cassette 
tape. 

The data was edited for gross errors or inconsistencies due to 
known instrumentation malfunctions. The criteria for retaining or 
discarding data from individual levels or for entire ten minute inter- 
val depended on performance check results obtained during measurements, 
obvious inconsistencies between levels and sequential times, and the 
availability of fluctuation data. After mean data were obtained for 
ten minute intervals and passed preliminary editing for obvious erro- 
neous values, it was processed for three applications, i.e., mean wind 
U, mean temperature T, and mean humidity qg. These values were plotted 
on 2-cycle logarithmic paper. Since U, T, and g are expected to vary 
logarithmically with height, a best fit straight line should describe 
the distribution of data points. Best fit lines were used to estimate 
gradient parameters and the Richardson number (equation 5a, Figures 10 
and 12). In general, the procedures were subjective and in some in- 
Stances, best fit lines with different slopes could be visualized for 


a given set of data points. Therefore, one criterion used was not to 
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give a single data point too much influence in determining the best fit 
line. Consequently, the line drawn represented a most probable posi- 
tion between data points without necessarily coinciding with any point 


(Figure 12). 


B. ANALYSES OF FLUCTUATION DATA 

Spectral analyses to obtain a and e values were performed on 
fluctuating velocity and temperature data obtained with single wires 
Oriented in the vertical normal to the mean wind. These data represent 
a point measurement and yield temporal or time descriptions of the 
fluctuations. Twenty-one minute segments of data originally recorded 
On magnetic tape were recorded into the EMR-Schlumber model 1510 digital 
Spectrum analyzer. Procedures for converting spectral values, obtained ~ 
with the analog spectral analyzer, to engineering units and for obtain- 
ing turbulence parameters from the spectra are described in the follow- 
ing paragraphs. 

A necessary procedure was to scale the spectral plots in order to 
relate rms input voltages to power spectral densities (PSD), variance 
per unit frequency. To obtain PSD levels, corresponding to rms voltage 
Input, calibrated scales charts had to be constructed. 

For purposes of the X-Y plot format of the analyzer output, the rms 
voltages were converted to y = 1094 (voltage) units and a graduated 
Scale was constructed enabling the logarithm of rms volts to be inter- 
polated from the spectral plots. The vertical scale, Y, values from 
the plots were adjusted for each spectrum as a function of input and 
Spectral gains. The values were then converted to PSD levels for use 
in calculating turbulence parameters. The relation used was 


(cal. level = Y.2 
S(f) {PSD units} = ae GHGAnaeni (10 ) (20) 
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Figure 13. Virtual potential temperature profiles. 


of 





where bandwidth is the frequency range divided by the number of 


channels 


Freq range 


BW = Number channels 


256 Hz/256 = 1 Hz 


The calibration level volts rms equals the voltage at Y=0. (1 volt rms 
for a 3.16 volt input setting). 

Amplitude scaling calibrations were accomplished using externally 
generated white noise. Signals with a 1 volt rms over a frequency 
range of 0.1 Hz to 1000 Hz, giving a PSD of 107° v°/Hz. Setting a 
3.16v (10 dB) input gain on the spectrum analyzer insures that 1 volt 
input corresponds to Y=0. An example of such a calibration plot is 
shown in Figure (15). 

The turbulence parameters - and e were obtained from the 
variance spectra on the basis of the universal formulae for the inertial 
subrange in wave number space. These expressions predict a -5/3 slope 
for variance spectra of both velocity and temperature when plotted in 
log-log format. Figures 14 and 15 are typical spectra considered in 
analyses. The velocity variance spectra have consistent -5/3 slope. 
However, temperature spectra often exhibited slopes slightly different 
from the expected value of -5/3. This feature of temperature spectra 
has been observed by others and the existence of "cold spikes" in the 
temperature traces has been postulated as a probable cause, Friehe 
(1976). 

Assuming -5/3 slopes for the variance specta, the intercept of the 
Slope with the 1 Hz frequency line was the power spectral density value 
used in computing the parameter of interest, e and one : 
38 
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The measured PSD value obtained from velocity spectra was converted 


to a spectral density in engineering units by the equation 


=) 
2. PSD 


S(f) oF 


(cm/s/volt)* , volt®/Hz = (cm/s)*/Hz 


(21) 


where C is the hot-wire calibration factor. Similarly, the equation 
used to convert PSD values obtained from temperature spectra to 


dimensional values was 


2 


by 


[Pop 


th 
ee” 
i} 


\e 


(°C/volt)~ « vol t*/Hz = °¢e 7H (22) 
where C, is the calibration factor for the temperature system. 

Since velocity and temperature fluctuations were measured at a fixed 
point in the flow, the resultant spectral estimates correspond to 
"temporal" frequencies, f, as indicated by the notation Sy {f) and 
S+(f) in the two preceding equations. To obtain ce and c,. values, 
temporal (f) and space (k) scales must be related in order to use 
equations (21) and (22). This is accomplished by Taylor's (1938) 


"frozen turbulence’ hypothesis, that is 
k = 2nf/U (23) 


where U is the mean wind at the measurement level. The term "frozen" 
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implies the turbulence pattern remains unchanged as it Sweeps past 
the sensor probe. 


The following identies are those used to relate temporal spectral 


estimates to spatial spectral estimates 


-c/3 p7e/3 


FS, (f) (24) 


nT} 
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Oe Ome (25) 


2 T 


where S(f) and S7(f) are spectral density values with units of 
(m/s)°/Hz and ce Hz respectively. The above equations are rewritten 
in the following forms to obtain e and ae from values of f, U, 


eet), and S-(f) 


U 
2/3 
c= (+ f+ s(t)’ (26) 
2/3 
i : i Me 6 SA (27) 
2 


where k = 2nf/U and C, = .5 and C, = .25 , empirically. 


C. ANALYSES OF PAIRED TEMPERATURE SENSOR 

e values were also estimated by using paired sensors for variance 
analyses. This analysis was based on the expression for the tempera- 
ture structure function, equation (2). 

The analysis iS one in which temperature differences between paired 
resistance wires separated by a distance r is subjected to variance 


analyses as indicated by the term IT" (x)-T' (xr) 1° On equation (2). 
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Voltages corresponding to the temperature differences, AT, as 
measured by the paired sensors were recorded on analog magnetic tape. 
During experimental periods, calibration voltages were ane recorded 
on the data tape for reconstruction of the signal. This calibration 
Signal was played back at the time of data analysis to provide a 
reliable check of the recording system performance in both record 
and playback modes. 

The variance or mean square of the temperature differences required 
to compute ae was obtained by analyzing the analog voltages with a 
Mean square voltieter (Thermo Systems Incorporated, model 1060) using 
a 30 second time constant for the averaging. The output of the mean 
Square voltmeter as well as a record of the original signal were recorded 
on standard strip charts. A strip chart section appears in Figure 17. 
The record of the original instantaneous signal provides an important 
means of spotting extraneous noises. The rms voltmeter-strip chart re- 
corder system was also calibrated using known input voltages. At least 
One input calibration voltage was larger than the largest expected signal 
on the data record and likewise one input calibration voltage was lower 
than the smallest expected signal. This procedure minimized any non- 
linearities in either the rms voltmeter of the strip chart recorder. 

Five minute averaged values of the rms voltages were abstracted from 
the strip chart and these values are utilized in a program for the IBM 
360 computer. This program converts the mean square voltages values to 
temperature difference variances values, from which ae was computed. 
A final check on the validity of the measured data was made by perform- 
ing spectral analyses on representative data samples from each experi- 


mental run. The spectrum analyzer, described in section B, was used to 
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obtain the spectra which were subsequently plotted and kept as a 

permanent record. Any harmonic noise appeared as a spike in an other- 
wise reasonably smooth curve. It should be emphasized that continuous 
calibration and data checking procedure were essential steps through- 


Out the analyses. 


D. HOT WIRE CALIBRATION 

In-situ calibrations of the velocity sensors were accomplished 
every 10 minutes during the experiment. In this procedure recordings 
were made of both the cup anemometer wind speed and the corresponding 


hot wire voltage output. The sensor wind speed is given by 


2 |i 


Vo=aU’° +b (28) 
where V is the hot-wire voltage output and U is the mean wind speed 
for any given level. The constants, a and b, are the calibration curve 
Slope and intercept respectively. The constants, a and b, are derived 
from the in-situ calibration curve. 


To convert the voltage PSD levels to velocity units requires a 


calibration factor given by 
i eae. (29) 


where C is the calibration factor in m/s /volt, v' is the voltage 
fluctuation, and U' the velocity fluctuation. Differentiating equa- 


tion (28) yields 


Ve (ayul/*a) v' (30) 
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which when substituted into equation (29) yields 
ba= avul/*/a (31) 


An example in-situ calibration curve is shown in Figure 18. The 


slope of the curve yields the value for a in equation (31). 
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A comparison of the observed overwater results and those obtained 
by Wyngaard, et al. (1971) appear in Figures 19, 20, respectively. 

In each case dots represent an individual data point. For the overwater 
results mean values over Ri intervals of 0.25 were determined, plotted 
in the center of the interval, and appear as small dots within circles. 
Error bars delineate one standard deviation from the mean within each 
interval. The number at the top of each error bar is the number of 

data points used in defining the mean in the interval. 

For the stable and unstable regime with 0.5 > Ri >-0.5 there is 
some correlation with the results of Wyngaard, et al. However, for 
Ri <-0.5 there is little agreement with the difference in results on 
the order of one order of magnitude. There is some agreement with the 
trend in the sense that it does not change significantly with 
instability. 

Many of the Richardson numbers were not within acceptable ranges, 
for positive, and negative Ri values. Results are not included for 
positive Ri values above 1.0 and negative values lower than -3.0. A 
predicted critical positive Ri value can be derived from equation (19), 
which gives a value of 0.2]. 

Scatter in the observed results could be attributed to the scatter 
in both the measured = and 30/dz values. Deviation of tempera- 
ture spectra from the -5/3 slope caused uncertainty in the ae esti- 
mates. The 00@/dz values are dependent on profile estimates whose de- 
termination was at times subjective as discussed in conjunction with 


Figure 21. 
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Figure 19. Wyngaard, et al. (1971) prediction’ curve. 
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Friction velocity values, U, , were calculated in two ways. In 
the first method (Method A) profiles of log (z) versus mean wind, U , 
were plotted to obtain gradient estimate. An example profile is 
Shown.in Figure 10. After the profiles were drawn, 5 and 10 meter 
level wind speeds were estimated from the profiles for each 10 minute 
period and these values were used to solve equation (14) for U, . 
These values appear in column A of Table I. 

Friction velocity values were also calculated from fluctuation 
data (Method B). First, equation (24) was solved to obtain a value 
for the turbulent kinetic energy dissipation rate, e« , at each level. 
Then the neutral prediction, equation (16), was used to obtaina U, 
at each level. Ideally U, at all levels would be equal for any 
given period since U, is assumed to be constant in the surface layer. 
The values, with about 10% exceptions, were within acceptable limits 
of each other. They were then averaged to obtain representative 
values for U, . Results from these computations appear in Table I 
column B. 

U, values were also calculated using wind speed dependent Cio 
values suggested by Cardone (1969) and a value of U_ taken from pro- 


files, at a height of 10 meters, using 


Us = UiprCi, (32) 


Cardone's Ci representation was for the marine boundary layer. 
These values are shown in Table I column C. 

The results are similar to those reached by Atkinson (19/6) 8 1 rer. 
Showing some disagreement between U, values obtained by the two 
methods, with the fluctuation data method more closely agreeing with 


Cardone Ss results. 
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The normalized dissipation rate KeZU, > in relation to Z/L was 
tested in the same manner as Garratt (1972). The results of Garratt 
(1972) and those obtained in this study are shown in Figures 22 and 23, 
respectively. The curve represents the function o,(2/L)-2/L esa lar 
Scatter was also observed in an investigation by Busch and Panofsky 
(1968). 

An examination was also made of the variation of the observed dissi- 
pation rate of turbulent kinetic energy, €« , with height for the differ- 
ent stability conditions. Results of this appear in Figures 24 and 25 
for stable and unstable conditions, respectively. The results from 
both stability regimes are compared to the -1 properties for neutral 
Stability. The results show that for all stability conditions a slope 


corresponding to near neutral conditions was observed. 
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Figure 23. Observed results. 
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VI. CONCLUSIONS 


Correlation of spectrally derived ae results with the stability 
Parameter, Ri, was poor for the data obtained in this study of the 
marine environment. Verification evaluations for friction velocity 
results from mean (profile) and turbulent (dissipation) data were satis- 
factory. Correlation of normalized e with the stability parameter 
(Z/L) showed the same scatter as results of others. The variation of 
€ with height exhibited results which support the near neutral pre- 
dictions for most cases. 

It is concluded that the turbulent field in the marine boundary 
layer is subjected to several anomalous effects which can cause the 
disruption of the inertial subrange. Thus, empirically, axpressione 
for describing overland transmission characteristics will undoubtedly 
have to be altered to provide operational determinations for utiliza- 


tion in the marine environment. 
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Period Level d6/dz es 3 Ri 
(M) °C/cm x 10 
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Period Level d6/dz T.3 Ri 
(M) 2C/em x 10 
1974 8-13-1742 4.19 -0.025 4.25 -0.45] 
6.60 -0.016 0.45 -0.711 
7.64 -0.013 0.63 -0.823 
lence -0.007 Oa -1.496 
9-18-1643 4.19 0.031] 0.29 0.308 
6.60 0.020 0.08 0.486 
7.64 0.017 - 0.500 
lomeg 0.009 Oa] 1.000 
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ees -Q.00] 0.96 -0.2/70 
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13.89 -0.001 Tag -0.330 
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VSa6¢ -0.003 Ogi -2.013 
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6.60 -0.003 Voll of OP ols) 
7.64 -0.003 - -0.180 
lSaee -0.001 1.14 -(). 328 
9-19-201 | 4.19 -0.010 0250 OFU 
6.60 -().006 0.8] 0.080 
7.64 -0.005 - 0.093 
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Period Level d6/dz C13 Ri 
(M) °C/cm x 10 
1974 9-19-2158 areal? -0.010 0.46 -0.202 
6.60 -0.006 UA -0.319 
7.64 -0,005 - -0.369 
13.89 -0.003 1.10 = ter || 
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13.89 — 20.054 1.08 -3.000 
9-20-0014 feakg -0.011 0.66 -0.461 
6.60 -0.007 (a5 | 1 =ORe/ 
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6.60 -0.002 os salle | Ss 
7.64 -0,002 - eile | Os 
13.89 -0.001 1.68 =Oneoe 
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